A geometrical approach to the dynamics of spinor condensates II: Collective modes 
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In this paper we study the linearized dynamics of spinor condensates using the spin-node formalism 
developed in Ref. [1]. We provide a general method to linearize the equations of motion based on the 
symmetry of the mean-field ground state using the local stereographic projection of the spin nodes. 
We also provide a simple construction to extract the collective modes from symmetry considerations 
alone akin to the analysis of vibrational excitations of polyatomic molecules. Finally, we will present 
a mapping between the spin-wave modes, and the wave functions of electrons in atoms, where the 
spherical symmetry is degraded by a crystal field. 



I. INTRODUCTION 

Spinor condensates have become a central theme in 
atomic physics since their initial creation [2-7]. Recent 
theoretical interests in spinor condensates have focused 
on topics such as dynamics near the insulating transition 
[8], metastable decay of currents [9], spin knots [10], and 
the anomalous Hall effect [11]. One particularly impor- 
tant aspect of spinor condensates is their free dynamics 
under a time-dependent Hamiltonian, about ground or 
metastable states. This aspect was the center of several 
experimental [6, 7, 12-15] and theoretical [16-22] stud- 
ies. These investigations, so far, were mostly confined to 
the simplest case of spin-one condensates. On the other 
hand, the wealth and intricacy of spinor condensates in- 
creases dramatically with increasing spins. For instance, 
the phase diagram of spin-two and spin-three condensates 
consists of four, and ten possible mean-field phases, re- 
spectively [23-28]. A feature which makes these systems 
even more interesting, is that the ground states exhibit 
a high degree of symmetry in its spin state, which is iso- 
morphic to lattice point groups [29-31]. In this paper 
we seek to utilize this symmetry in the study of the free 
dynamics of spinor condensates. 

In a companion paper [1], we utilized the spin- node ge- 
ometrical description of the ground state of spinor con- 
densates and generalized it to study the hydrodynamics 
of these systems. Our goal was to introduce hydrody- 
namic equations of motion that preserved the geometri- 
cal structure of the spinor degrees of freedom, in a way 
that the Gross-Pitaevskii equation (GPE) cannot - the 
GPE requires the choice of a fixed axis of spin quanti- 
zation, and it therefore hides the symmetries of the spin 
wave function. In Ref. [1], we expressed hydrodynamic 
equations of motion using the density, superfluid veloc- 
ity, and the spin-nodes as our basic degrees of freedom. 
In addition to the Euler equations, which describe mass, 
momentum, and energy conservation, we obtained 2F 
Landau-Lifshitz equations for the dynamics of the spin- 
nodes. Furthermore, our derivation gave a natural gen- 
eralization of the Mermin-Ho relation which connects the 
vorticity in a ferromagnetic spinor condensate with the 
Pontryagin index of the order parameter. We were able to 



give explicit expressions for all dynamical quantities for 
spin-half and spin-one condensates. However, we found 
that the full expressions become progressively more com- 
plicated with increasing spin, and become impractical for 
large F. 

In this paper, we study the dynamics of spinor conden- 
sates for general values of F, by looking at the low energy 
properties of the equations of motion obtained in Ref. 1. 
Our goal is to derive the spin-node description of the low 
lying spin-wave excitations near the mean-field ground 
state. In order to do this, we derive the linearized equa- 
tions of motion for the spin-node locations, which allows 
us to extract the small oscillation spectrum from symme- 
try alone, in a fashion resembling the vibrational-mode 
calculation for polyatomic molecules [32, 33]. Using this 
method we are able to give simple expressions for the vi- 
bration eigenmodes and energy spectrum. In addition, 
we derive a correspondence between the low lying excita- 
tions of the spinor-condensates, and atomic orbitals sub- 
ject to rotational symmetry due to crystal-fields, which 
reflect the symmetry of the spinor-condensate ground 
state. 

The paper is organized as follows. Sec. II provides a 
brief summary of the main results from Ref. [1] which will 
be used in this paper. In Sec. Ill we derive the linearized 
equations of motion about the mean-field configuration 
in terms of the spin-node formalism. Finally, in Sec. IV 
we demonstrate how to use symmetry arguments to com- 
pute the spin- wave excitations, and give a prescription to 
obtain closed form expressions for both eigenmodes and 
cigenencrgies of the low-lying spin- waves. 



II. SPIN NODES AND HYDRODYNAMIC 
EQUATIONS 

In this section, we summarize the notation and results 
from Ref. [1] which will be necessary for the development 
in this paper. We first discuss a single F = g spinor 
and its associated vector triad, and then describe how 
the mean-field ground state of a spin-i* 1 condensate is 
described in terms of 2F spin nodes. Finally, we give 
a brief summary of the hydrodynamic equations for a 
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spinor condensate with general spin F. 



B. Spin node representation for spin-F condensate 



A. Spin- half spinors and vector triads 



i spinor can be described in terms of a 



A single F 

solid angle il = [9, 0), corresponding to the direction on 
the unit sphere in which the spin is pointing: 



To describe a spin-i 7 condensate, we begin by sepa- 
rating the wave function into a piece corresponding to 
the overall density and phase, and a piece describing the 
local spin state. We write 



|0) 



i<j>/2 



IT) 



,-i<j>/2 



(1) 



This specifies the wave function uniquely, up to an overall 
phase. In what follows, we adhere to the phase conven- 
tion above, which corresponds to a particular choice of 
gauge. However, the equations of motion that we will 
derive will be shown to be explicitly gauge invariant. 
The time-reversed partner of |f2) is, 



"0a = IpXa 

where Xa is a normalized spin-i 7 spinor 

^X* a Xa = (X\X) = 1 



and the superfluid density is 



H 5 



(9) 



(10) 



(11) 



sin 



IT) 



e 4 ^ 2 |I) . (2) 



which satisfies (f2*|f2) = 0. In addition, we define 

\n x ) = -L(\n) + \n t )) (3) 
\n y ) = -L(|n) + i|n t ». (4) 



We now turn to the spinor wave function \a = ( a lx)i 
which is the main focus of this paper. 

A general (non-normalized) spin-i 7, spinor \ can be 
represented as a totally symmetrized product of 2F spin- 
half spinors |Qi) . . . \{1%f) as i 



in) = |Oi 




(12) 



n = 2 <0| P |0) , 
e x = 2{fL x \F\n x ) 
e y = 2(Q y \F\Q y ) 



(5) 



where the sum is over all permutations a over 2F indices 
(note that the spinor (12) is not normalized to unity, 
These states allow us to introduce an orthonormal triad anc j we rese rve the notation \ X ) = , |n> for normalized 
of unit vectors (n, e x , e y ), defined by _ _ v< n ^> 

spinors). Thus, we can think of a general spin-i 7 spinor as 

a collection of 2F indistinguishable spin-half spinors. In 
what follows, we call these spin-half spinors 'spin nodes'. 
Thus, a spin-i 7 spinor is built out of 2F spin nodes, each 
one of which has an associated vector triad, as discussed 
in Sec. II A. The relationship between spin nodes and 
the reciprocal vectors of Ref. [29] is discussed in detail in 
Rcf. [1]. 

There are a few special spinors and that will be espe- 
cially useful in what follows. First, a spin-coherent state, 
denoted |(f2) 2F ), is a spinor in which all spin-nodes are 
chosen to be the same (i.e. they all point in the same 
direction), 



Here, n is simply the unit vector along which the spinor 
points, whereas e x and e y are the two directions perpen- 
dicular to it. There is an ambiguity in the definition of e x 
and e. y corresponding to the gauge choice in the defini- 
tion of |f2). Gauge invariant quantities can be expressed 
in terms of n only, without reference to e x and e y . 

In Ref. [1] we derive a number of useful identities for 
the spin-half spinor and its associated triad. For instance: 



and 



a a = i (D,\d a Cl) 



(f2'|9 a f2) = — e+ • <9 a n. 



(6) 
(7) 

(8) 



(n) 2 ") = \n...n) 



(13) 



We next define the spin state corresponding to | n) with 
its ith component time-reversed. We denote these by 



(14) 



In the second expression, we have introduced the complex 
vectors e± = ± ie y . In the first equation, we use the 
notation a a to reflect the role of the quantity i (£l\d a Q) as 
the vector potential associated with the gauge symmetry 
of the spinor. 



Finally, we define the projection operator V to be 

P = 1-|X><X|. (15) 

See Ref. [1] for a more thorough discussion of the prop- 
erties of the spin node representation, including a careful 
treatment based on the Schwingcr boson approach. 
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C. Hydrodynamics for general spin-F condensates. 

In Ref. [1] we derived the hydro-dynamic equations of 
motion for general spin-F 1 . The first two equations of mo- 
tion, the mass continuity equation and the Eulcr equation 
are 



d t p = -V • (pv) 



and 



D t w = e + (v x b) - V 



2V5n 



(16) 



(17) 



p 2 2^p 

The effective electric and magnetic fields follow from the 
field tensor f a( g constructed from a a = i (x\d a X)- 

To obtain the Landau-Lifshitz equations, we contract 

the GPE with ( (n*) . Doing this gives 



-d Q log 



4(^|v 2 Oi>-(o*|a a Jii)^ 



a* n n 



(n*) 2F \(n\v int \n) 



R'l%) 



in 



(18) 

where the subscripts of the bra's and ket's in the last term 
denote how the inner product is to be evaluated: ket 1 
(2) is contracted with bra 1 (2), and signify the state of 
one of two interacting particles. In this expression we 
have introduced the quantities A^, 

\ i = (2F)\H(n j \ty). (19) 

See Ref. [1] for a more complete derivation of this ex- 
pression. 

While the first term in Eq. (18), 



d f n. 



is the inertial term for the spin- node n^, the right hand 
side, and the last term of Eq. (18) in particular, should 
serve the role of torques, projected onto ej + . As we will 
show in the next section, the matrix element of Vint is in- 
deed related to a derivative with respect to the spin-node 
coordinates of a potential energy function. Specifically: 

P 



a* (nin) 



, 2F. 



jtfp <n|n) 



;ni Vi„t. in) .in ). = 



[e J+ -V nj V({ ni })] 



(20) 



where V({rij}) = (Vint) is the expectation value of the 
energy of a spin configuration with spin nodes {n^}, and 
AT. is a matrix which projects the torques due to spin- 
node j and the motion of spin- node i. The matrix A 
and its inverse are defined below in Eq. (29). Eq. (18) 
provides a natural starting point in the analysis of the 
linearized equations of motion which will be developed 
in the following section. 



III. LINEARIZED EQUATIONS OF MOTION 
FOR ARBITRARY SPIN-F CONDENSATES 

As suggested from the equations of motion of the spin- 
one and higher condensates in Ref. 1, the geometric repre- 
sentation of the equations of motion yield rather compli- 
cated results. Nevertheless, this formalism regains its ap- 
peal when linearized about particular mean-field ground 
states. Then the hidden point symmetries of the ground 
state become apparent, and can be used to describe the 
linearized dynamics of a condensate. Below we derive 
the small oscillation description of general spinor con- 
densates. 



A. Linearized equations of motion from the GPE 

Parting ways from the attempt at a general descrip- 
tion of spinor condensate dynamics, we now turn to the 
vicinity of a uniform mean-field ground state. For the 
ensuing discussion, we will denote quantities to be eval- 
uated in the mean-field ground state with overhead bars. 
For instance, the density can be written by expanding 
about the mean-field state as 



p = p + 5 p. 



(21) 



We will first concentrate on the equations describing the 
density excitations. To linearize the equations of motion, 
derived in Sec. II C, we can drop terms which involve 
derivatives acting on two different quantities. Doing so 
leads to the following two equations describing the den- 
sity fluctuations: 



and 



d t p = —p V ■ v, 



2VW V 2 p 

d t v= -V —p+ — 

P 2p 



(22) 



(23) 



Note that terms describing the spin degrees of freedom 
(e.g., the effective electric and magnetic fields) have com- 
pletely dropped out of these equations from lineariza- 
tion. Computing the excitations from these equations is 
straightforward and gives the familiar Bogoliubov mode 
describing density fluctuations. 

Let us now focus our attention on linearizing the 
Landau-Lifshitz equations for general spin written in 
Eq. (18). Since the process of linearization separates the 
equations for spin and density fluctuations, to simplify 
the notation in what follows, we will scale the density of 
the uniform state to one, po — > 1. When linearized, the 
Landau-Lifshitz equations for general spin become 



= — (ft 

2 x 



4| V 2 ^) 



(24) 



A*(n|n) iX 



m 



2F. 



'niVintin) in 
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In the above, as before, we have used overhead bars to 
denote quantities evaluated at their mean-field configu- 
ration. 

To understand the dynamics of Eq. (24) it is useful 
to introduce variables to describe small deviations of the 
spin nodes from their mean-field values. To this end, by 
using the identities established in Sec. II A, we introduce 
the set of 2F complex variables {z^ 



1_ 
2 



1_ 



6m, (25) 



where = ri; + Siii . Note that in the mean- field states 
we have Zi = for each spin node since the vectors ei + 
and hi are orthogonal. This set of variables can be seen 
to be the local stercographic projection of n, onto the 
complex plane for small displacements, and will be very 
useful in the following analysis. Moreover, in our gauge 
convention, z% is given in terms of displacements along 
the zenith and azimuthal directions from the spherical 
coordinate system: 



B. Perturbative expansion of the spin interaction 

An essential element in the behavior of spinor conden- 
sates is the spin interaction term VJ nt . It is the minimiza- 
tion of this term that yields the mean-field ground states, 
and its curvature that determines the normal excitations. 
These curvatures can be easily and directly extracted in 
terms of specific matrix elements, as we show below. 

To expand the spin interaction energy about a mean- 
field ground state (denoted with an overhead bar) we first 
need to understand how to perturb a spinor about a fixed 
value. The following spin-half identity proves to be quite 
helpful: 



\sn) = |o) (n\5Q) + \n t ) 



(31) 



where we used the resolution of the identity in terms 
of fl) and its time reversed partner, and the definition 
of z = (Qt \SO) as in Eq. (25). Now, if we apply the 
variation to a general spin-F spinor |fi) = fi) + 6 |fl), 
we obtain to linear order 



Sn ■ 9 + iSn ■ (p. 



(26) 



2F 



■IF 



Using these variables, the linearized Landau-Lifshitz 
equations become 



id t Zi = -\v 2 z l + 



,2F. 



[silvan) in 



A*(n|n) 



(27) 



The kinetic pieces in the GP equations are most natu- 
rally described in terms of the original spinor wave func- 
tion, ijj a , and are not simplified by the symmetry of the 
mean-field ground states. Nevertheless, Eq. (27) demon- 
strates that near mean-field ground states the kinetic 
terms still acquire a simple form. Interestingly, the ki- 
netic parts in the spin equations of motion, (27), do not 
disclose the fact that the variables {z^ describe spin- half 
components of a spin-F state. This fact is reflected only 
in the spin interaction term. In the following section, 
we will see that this spin interaction can be expressed in 
terms of a derivative with respect to the z* variables. In 
particular, the equations of motion will be shown to be 



1 



idtZi = --v 2 z 4 + (nln)^Ir. 1 
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(28) 



where 



A7} ee 



2F 



(29) 



Thus, the spin interaction derives from a sum over 
"torques," 
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(30) 



s |n) = |n) + \ Tifl ) ■ 



(32) 



where |T^n) is |n) with its ith entry time reversed (see 
Sec. II B). Since /f2j|<5fij) is imaginary, the first term, 
which docs not directly depend on z, must drop off when 
considering the variations of real quantities. For instance, 
the first order variation of the normalization is: 



IF 

(n|n) = Y {(n\Tin) z, + (i;n|n) 



Using Eqns. (32,33) one finds to lowest order 



d in) 



where 



_ V\TM) 

dzi (n|n) _ (n|n) ' 



|n) (n| 



(n|n) 



(33) 



(34) 



(35) 



Such an expression is useful in evaluating derivatives 
of the spin interaction energy as in Eq. (28). In gen- 
eral, derivatives with respect to z* will act on bras while 
derivatives with respect to Zi will act on kets. 

We will now establish the equivalence between 
Eqns. (27) and (28). One can use Eq. (34) to evaluate 
the derivative of the interaction energy 



_d_ 

dz~" 



V in t = 



rr^n^njKrtjnjjn 
(n|n) 2 



(36) 



which is correct to linear order. The subscripts of the 
bra's and ket's denote how the inner product is to be 
evaluated: ket 1 (2) is contracted with bra 1 (2), and 
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signify the state of one of two interacting particles; sim- 
ilarly, the projection V\ operates only on the degrees of 
freedom pertaining to particle '1'. Then using the ex- 
pression for A^ 1 and the relation (derived in Appendix 
A) 



E 



mf F )(T t n\ 



(37) 



one immediately finds for the last term in Eq. (28) 



(41), only one projector in needed in Eq. (42), so Ti can 
be omitted. 

While these results for the spin interaction seem in- 
volved, they are directly expressed in terms of easily- 
constructed matrix elements evaluated at the mean-field 
ground state. Furthermore, these matrix elements obey 
the point symmetry of the ground state at hand, and 
thus have stringent constraints. Eq. (42) therefore pro- 
vides us with direct expressions for the matrix elements 
appearing in the linear spin- wave expansion of the spinor 
condensate. 



•14, 



2F, 



:m Vi„t in) jn 



a; <n|n) 



which is the last term in Eq. (27). 



(38) 



1. Second order expansion of the interaction energy 

Since we are interested in small oscillations about equi- 
librium, we would like to express the interaction energy 
expanded about the mean-field state to quadratic order 
in the z variables. This can be formally written as 



Vint — Vint 



d 2 V m 



2 ^ dzidz; 



dz*d Zj 



d 2 Vi 



int * * 

Zt Za 



dzfdz* 1 3 



(40) 



where the terms involving derivatives of Vi„t are to be 
evaluated at the mean-field ground state. We can now use 
Eq. (34) to evaluate these derivatives of the interaction 
energy. Note that terms where two derivatives act on the 
same bra or ket will vanish since 



p la <n| Vmt =o, 



(41) 



which happens since Ti = at the minimum of the spin 
interaction, so that (n|Vj n t|0 ) |S7 ) oc \ti ) . We then 
readily obtain the following quadratic form for the spin 
interaction energy (dropping the Vj n t term): 



int 




OI^IVint^l^O^llTj-n) 

2(n\n) 2 
T l n\P l2 (n\v int \ti)P 1 \T J n) i 
(n|n) 2 

i (T z n\v l2 (T 1 n\r 2 v int \n) 2 \n) i 
2(n\£i) 2 



z iZj 



(42) 



Here is the projection operator which only acts on 
states denoted with subscripts 1 or 2 respectively. While 
the form above is written symmetrically, following Eq. 



C. The Lagrangian of spinor condensates near 
equilibrium 

The equations of motion can be arrived at by expand- 
ing the spinor condensate Lagrangian to quadratic or- 
der in the z variables, and computing the corresponding 
Euler-Lagrange equations. As we saw before, to this or- 
der, the density excitations decouple from the spin ex- 
citations. Thus, to simplify the analysis, we will fix the 
density and scale it to one, and work in the incompress- 
ible regime. The Lagrangian for a spin-_F condensate in 
the incompressible regime is 



C = a t 



1 



(V0 



1. 



:T - V m 



(43) 



where Vi„t is the spin interaction potential. In expanding 
this Lagrangian to second order, we first consider the spin 
Berry's phase contribution 

Note that the kets and bras involving time derivatives 
are necessarily first order in variation from the mean-field 
state. Thus we consider the following quantity expanded 
to first order about the ground state 



|n) 



(fiin) 



\sn) 
(n|n> 

V\Sil) 

(n|n) 



(n|n) 
(sn\Ct) 



?5((n\n)) (45) 



(46) 



Inserting this into the expression for the spin Berry's 
phase (44), and dropping terms that can be written as 
total time derivatives (which do not contribute to the 
dynamics) one finds 



n, 



(n\r\d t n) 



(47) 



We can then insert into Eq. (47) the expressions for 
the expansion of |f2) to linear order in the z variables 
given in Eq. (32) to directly obtain 



(n|n) 



ij 



AndtZj 



(48) 
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where 



Aij = (TiSllPlTjSl) 



(49) 



which is the sought-after relation. The proof that A de- 
fined here is in fact the inverse of the expression given in 
Eq. (29) is given in Appendix A. The hermitian matrix 
A gives the canonical commutation relations between the 
z variables. To directly compute the matrix elements of 
A is cumbersome because each involves a Wick expan- 
sion of (2F)l terms. On the other hand the expression 
for A^ 1 given in Eq. (29) is readily computed since it in- 
volves evaluating overlaps between spin-coherent states. 
Thus, in practice, to construct the matrix A it is easiest 
to first construct A -1 and then compute its inverse. 

Proceeding along very similar lines as above, one can 
expand T to second order in the z's. One finds 



{d aX \V\d aX ) 



1 



(n|n) 



^daZ-AijdaZj. (50) 



Finally, we note that the term involving the superfiuid 
velocity v = V# — a in the Lagrangian will not contribute 
to the linearized equations of motion. We are now in a 
position to vary the Lagrangian Eq. (43) as a function of 
the z's to find the linearized equations of motion. These 
read 



iAijdfZj 



l -A l3 ^z ] + {Sl\Sl)^ (51) 



(repeated indices are summed over) . It is straightforward 
to see that this is the same as Eq. (28) which was obtained 
directly from linearizing the GPE contracted with time- 
reversed coherent states. 

Since A is a hermitian matrix, it is diagonalizcd by a 
unitary transformation 



A = UAUK 



(52) 



where A is the diagonal matrix consisting of the eigen- 
values of A. It is therefore convenient to define a new set 
of ^-coordinates as 



(53) 



Note that in terms of these coordinates, the Berry's phase 
assumes a simple diagonal form 



(«|n) 



(54) 



Furthermore, the equations of motion have the simple 
form in these coordinates: 



id. 



i_ o (nln) dv 
' i = l v ^TT^ (55) 



This has the form of a time-dependent Schrodingcr equa- 
tion for the Wi parameters. 



IV. NORMAL EIGENMODES, SYMMETRY, 
AND GROUP THEORY 

The most appealing application of the linearized equa- 
tions of motion developed in the previous section is to ob- 
tain the normal excitation modes and energies of spinor 
condensates having a hidden ground state symmetry. As 
we show, it is nearly sufficient to diagonalize the matrix A 
[defined in Eq. (29)] in order to obtain the eigenmodes of 
the spinor-condensate. This can be done solely by using 
the symmetry of the hidden symmetry of the mean-field 
state. 

Below we first demonstrate the use of the linearized 
equations of motion on the cyclic state without fully uti- 
lizing the symmetry in Sec. IV A, and obtain all eigen- 
modes and eigenfrequencies using the variables defined 
in Sec. IIIB1. Next, in Sec. IV B, we demonstrate 
how from the point group of the hidden symmetry of 
the mean-field ground states, we can compute the nor- 
mal modes alone (but not energies), using the example 
of the spin-three state where the spin-nodes are arranged 
at the vertices of a hexagon. Finally, in Sec. IV C, we 
show how to directly construct the vibrational and ro- 
tational eigenmodes from spherical harmonics, by con- 
necting the problem at hand to that of degeneracy lifting 
of electronic atomic orbitals. This method circumvents 
the arduous group-theory foot work, by using the well- 
known properties of atomic orbitals under crystal fields 
that break rotational invariance. 

The general motivation of the discussion below is that 
group theory analysis can be applied to obtain the nor- 
mal modes in spinor condensates much like the analy- 
sis of the vibrational frequencies of polyatomic molecules 
[32, 33]. The "atoms" (or spin nodes) in our case, how- 
ever, are confined to the surface of the unit sphere, and 
the displacement of each spin node is a two-dimensional 
vector (parameterized by the real and imaginary parts of 
the z variables). The first step in a symmetry analysis 
is to construct the transformation rules of the 2-d dis- 
placement vectors under point group symmetries. These 
transformation rules are a reducible representation of the 
symmetry group, and can then be broken down into its ir- 
reducible representations (irreps). The modes that trans- 
form according to the irreps are the eigenmodes of the 
system. 

Before we begin the analysis, a note on mode multi- 
plicity is in order. Naively, one might expect that the 
procedure in the previous paragraph will give (2F) x 2 
normal modes due to the two basis vectors per spin node. 
This situation would arise if the transformations we con- 
struct transform the 2F x 2 real coordinates, and are 
therefore 4F large real reducible representations of the 
symmetry group, resulting in 4F modes. While this is 
the case for real atoms, where the displacement vectors 
are also associated with conjugate momenta, the spin- 
nodes displacements do not have independent conjugate 
momenta. From Eq. (43) and (44) we see that the com- 
plex displacement zi is actually canonically conjugate to 
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t)z, 



the two-dimensional displace- 



ments arc both the coordinate and conjugate momenta, 
and hence there are only 2F eigenmodes in a spinor con- 
densate. Qualitatively, this is a situation reminiscent of 
a massless particle in a magnetic field, where the x and 
y coordinates are canonically conjugate coordinate and 
momentum. Indeed, constructing real 4F dimensional 
representations of the symmetry would result in two du- 
plicates of the spinor-condensate's eigenmodes. This du- 
plicity will become evident when the eigenmodes are writ- 
ten in terms of the complex Zi's: half the normal modes 
will differ from the other half through a complex multi- 
plicative coefficient. 



A. Spin- two cyclic state 

As our first example, we consider the cyclic state which 
is a possible mean-field ground state having the symme- 
try of a tetrahedron for the spin-two problem. We will 
expand the interaction energy to quadratic order about 
this mean-field ground state to compute the energies of 
the normal excitations. The spin-two interaction energy 
can be written in the simple form [25, 26] 



1 



1 



<*m d + -0\ ( X t\x) I' 



(56) 



where a and (3 are functions of the scattering lengths, 
and 



m = (x|F|x) 



(57) 



For the mean-field cyclic state, this spin interaction en- 
ergy conveniently vanishes VS„t = 0. In the following we 
will expand this energy to quadratic order. 

We first construct the symmetry matrix A for the cyclic 
state. We take the orientation where the spin nodes are 
at (in cartesian coordinates) 



1 



ni 



n 3 



n 2 



-1,1) 



(1,-1,-1), 



n 4 



v/3 



(-1,1,-1). 



(58) 
(59) 



With the spin-half spinors corresponding to these spin 
nodes the matrix A -1 can be directly constructed using 
the expression involving overlaps of time-reversed coher- 
ent states in Eq. (29). Using our gauge convention, this 
is found to be 




A~ 



This then can be inverted to obtain 



A=- 



(60) 




(61) 



Recall that directly constructing the A matrix is cumber- 
some since its elements involve Wick expansions having 
(2F)! terms. The eigenvalues of this matrix are found to 
be Eig(A) = (A~i, A 2 , A 3 , A 4 ) = (8,8,8,^). This matrix 
can be written in a revealing form as 



A = 81 W4U4 



(62) 



where U4 = 1, — 1,— 1) T is the eigenvector of A cor- 
responding to eigenvalue A4 and / is the identity matrix. 
An eigenmode will necessarily diagonalize the A matrix 
as well as the entire equations of motion, and therefore 
we already gleaned one eigenmode: U4, which will turn 
out to be the optical mode. 

The three modes orthogonal to U4 are associated with 
SO(3) rotations. With this in mind, we construct these 
three eigenmodes as the vectors arising from infinitesimal 
rotations of about the cartesian axes, x a . A rotation 
by angle 8rj about the x a axis produces the following z^s: 

Zi(Sr]) = 5r] (± a x n.A • e i+ = iSij e i+ ■ x Q . (63) 

Thus the eigenvectors u a are: 



V 1 

1 /; 



(64) 



It is now clear how to write the transformation into the 
cigen-coordinatcs defined generally in Eq. (53): 



(65) 



Due to the high symmetry of the tetrahedron, these mode 
are also degenerate. In general, the set of coordinate 
vectors x a should be taken to be the principal axes of 
the mean-field configuration. 

Next we use this matrix to expand the interaction en- 
ergy. We first consider the linear order variation of the 
spin moment m. Note that since m = in the ground 
state we have (ft\ F \fl) = FV \fl). Then by insert- 
ing the identity for V given in Eq. (37) one finds 



<5m = 



1 



2(n|n) 



]T(e 4 _A 



(66) 



We now write the vectors ei + in the basis of unit vectors 
along the three cartesian coordinates 



^2(e l+ ■ x Q )x a 



(67) 



which immediately reduces them to the complex conju- 
gate of the degenerate eigenvectors u a (with eigenvalue 
A = 8). The fact that all of the eigenvalues are the same 
is due to the high symmetry of the tetrahedral state. 
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TABLE I: The character table ol the group D§h using the 
notation of [32]. The last row gives the characters of the 
reducible representation F constructed from transforming the 
displacement vectors of the hexagon (see text). 



FIG. 1: Normal modes of the cyclic state. Mode (a) is the op- 
tical mode corresponding to pure displacements in 104. Modes 
(b), (c), and (d) are gapless modes corresponding rotating 
about the x, y, or z axes respectively. The axes of rotation 
for these modes is shown. 

With this basis one finds for the expansion of magneti- 
zation the simple expression 

\ I / i ol—1 
3 

= a/6^x«K+0 (68) 

where we have expressed the final result in terms of the 
w- variables (defined in Eq. (65). In deriving the above ex- 
pression, we have explicitly used the values for the eigen- 
values of the A matrix and the normalization constant 
^J7|f2) = |. The three parameters of Wi occuring in 
Eq. (68) correspond to rotations about the three carte- 
sian axes as shown in Fig. 1. 

Similar analysis can be performed on the second term 
in the spin interaction for the cyclic state. Without show- 
ing the details, it is found that 

5( Xt \x) = 2V2w*. (69) 

With these expressions we can now write down the spin 
interaction energy expanded to quadratic order which 
reads 

3 

(0|n) V s = Mm +w*) 2 +2/3A 4 K| 2 - (70) 

With this expansion of the interaction, Eq. (55) can be 
directly used to compute the energy of the normal excita- 
tions. Four Bogoliubov modes (note we are neglecting the 



density mode) are readily obtained. One finds three gap- 
less spin waves of dispersion Ef. = y^EkiEk + 4a) in addi- 
tion to an optical mode having dispersion E^ p = Ek + 2(3 
(where Ek is the free particle dispersion). 

Quite generally, the eigenvectors of the matrix A yield 
the displacements of the z variables corresponding to 
each of the eigenmodes (see, e.g., Eq. (63). In case 
of degeneracy, it is the interaction terms, discussed in 
Sec. IIIB1, that determine the correct diagonalization 
of the degenerate subspace in the matrix A. In the case of 
the cyclic state, the first three modes have displacements 
that correspond to rotations about three orthogonal axes. 
The final mode z oc u 4 corresponds to the optical excita- 
tion discussed above, and its displacements are depicted 
in Fig. 1. This procedure simplifies the standard Bo- 
goliubov method [34] considerably; we extract the eigen- 
modes solely from the A matrix, which, as we show next, 
can be obtained from symmetry considerations. 



B. Spin-three hexagonal state 

Let us now describe how to obtain the normal modes of 
a spinor condensate by using symmetry arguments alone 
in a more complicated setting. Once having the eigen- 
modes, however, we must note that to obtain the ener- 
getics and dispersions of these modes, analysis of the mi- 
croscopic Hamiltonian is still required. Our analysis uses 
group theoretical arguments similar to those used to de- 
termine the vibrational modes of polyatomic molecules 
[32, 33]. We illustrate the method through the nontrivial 
example of the spin-three state having the symmetry of 
the hexagon, which is a candidate for the ground state of 
52 Cr condensates [27, 28]. 

The hexagon belongs to the point symmetry group D§h 
whose character table is given in Table I. In this table we 
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FIG. 2: Normal modes for the hexagonal configuration of 
the spin-three condensate. Vectors moving into and out of 
the plane are denoted with "-" and "+" respectively. By 
multiplying the set of parameters {zt} corresponding to these 
displacements by a factor of i, one can identify v\ = V2, vz — 
v&, V5 = V7, v<i = vs, v§ = v\\, and fio = v\%. The modes 
v\ — V2, V5 = V7, and v§ = v%, correspond to Goldstone 
excitations due to the broken spin symmetry, while all other 
modes are optical and gapped. The mode fio = u\2 has a set 
of displacement vectors with lengths differing by a factor of 
two. 

use the standardized notation for the symmetry operators 
and irreducible representations [32]. To every spin node, 
we attach two displacement vectors parameterized by the 
real and imaginary parts of the Zj's introduced previously. 
Such displacement vectors are always parallel to the sur- 
face of the sphere. We can construct matrices which 
describe how this set of 2 ■ 2F = AF vectors transform 
under each of the symmetry operations. It is easy to then 
see that this set of matrices T = {Mi} form a (reducible) 
representation of the symmetry group. While these large 
AF x AF matrices are cumbersome to write down, their 
characters (traces) can be obtained by inspection. For in- 
stance, only spin nodes which are mapped to themselves 
by a particular symmetry operation will contribute to the 
character of the matrix describing this symmetry opera- 
tion. The last row of Tab. I gives the characters of each 
of the matrices Mi forming T. 

One can then invert the character matrix given in 
Tab. I to see how Y can be decomposed into combina- 
tions of irreducible representations. The result is 

r = A 2g + B 2g + E lg + E 2g + A 2u + B 2u + E lu + E 2u . (71) 

In the typical notation [32] ^4's and B's denote one- 
dimensional irreducible representations while £J's denote 
two-dimensional irreducible representations. The normal 
modes form the basis of each of these irreducible repre- 
sentations [32]. For two-dimensional irreducible repre- 
sentations, there is some ambiguity in picking the two 



basis functions. For simplicity, we picked the particular 
displacements which are all in-plane or all out-of-plane to 
form such bases. The AF modes corresponding to each 
of these representations is given in Fig. 2. As usual, the 
modes corresponding to two-dimensional representations 
are degenerate. 

Once we know the irreducible representations involved, 
we follow standard group theory, and construct projec- 
tion operators for the modes in these irreducible repre- 
sentations. A general displacement of the spin nodes Q 
can be decomposed into a superposition of modes forming 
bases for each irreducible representation as 

Q = £>(r i )Q (72) 

i 

where the operator V(Ti) projects into the irreducible 
representation Fj. Such projection operators can be writ- 
ten explicitly as 

^( r *) = tEx (ri) (.9) D (9)- (73) 
g 

Here, X {9) i s the character for the irredicuble rep- 
resentation Ti corresponding to group element g, ii is 
the dimension of the ith irreducible representation, and 
h are the number of elements in the symmetry group; 
D(g) is the representation of group element g in the spin- 
nodes displacement basis. For the hexagonal state of the 
spin-three condensate, this projection confirms the eigen- 
modes depicted in Fig. 2. 

As mentioned above, unlike molecular normal modes 
where the atoms oscillate linearly about the equilibrium 
positions, the spin nodes will rotate along ellipses about 
the equilibrium configuration. This allows us to cut the 
number of modes given Fig. 2 in half. Specifically, by 
multiplying the displacements {z{\ by the phase factor 
of i, we identify v\ = v 2 , vj, = v±, v§ = 1/7, v§ = v%, 
v<z = fix, and j^io = ^12- Because of rotational invariancc, 
the aspect ratio of the ellipses for the three spin rotational 
Goldstone modes v\ = v 2 , v§ = V7, and v§ = is s will 
be zero. Finally, we identify the three remaining modes 
v?, = U4, v% = Vn, and v\q = v\ 2 with gapped optical 
modes of the hexagonal spin-three condensate. 

Thus, for this spin-three problem, by symmetry argu- 
ments alone we have identified the IF = 6 spin modes 
(three of which are Goldstone modes). These modes 
along with density mode give the complete spectrum of 
normal modes for the spin-three hexagonal condensate. 

C. Connection to atomic orbital theory and 
spherical harmonics 

The treatment above makes the construction of low- 
energy eigenmodes of spinor condensates geometrically 
intuitive, and illustrates how to directly use the machin- 
ery of group theory. In addition, however, it is possible 
to make use of the close relationship of the symmetry 
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group and the underlying full 5*0(3) rotational symme- 
try (such a connection was explored in the context of 
equilibrium spinor-condensates in Rcf. [30]). Once this 
connection is made, we will be able to simply map the 
already well-developed theory of crystal-field splittings of 
atomic orbitals to the problem of eigenmodes of spinor 
condensates. 

The connection between the z-represcntation of small 
oscillations as in Sec. IV A and spherical harmonics can 
be deduced from the transformation rules of the vector z% 
under the relevant point-group. On the one hand, a sym- 
metry operator in the z-representation, Dij(g), will per- 
mute the entries Zj, as the symmetry operation g would 
the spin-nodes. On the other hand, each Zj is a two- 
dimensional vector written in terms of a complex num- 
ber with respect to a particular basis pair, e x , e v , which 
are functions of the location of the spin-node n on the 
unit sphere. Therefore, the operator Dij(g) also contains 
phase factors, e lXij ( g \ which serve to rotate the basis vec- 
tors. So, in general, the structure of symmetry operators 
in the z-basis is 



(74) 



1 2F) 

where A\j (g) is an element of the 2F permutation 
group corresponding to a rotational symmetry of the spin 
nodes. 

By exploiting the above transformation structure, we 
can systematically construct bases of the symmetry 
group over C 2F from the bases of rotational symme- 
try, namely, spherical harmonics, Yi m {9, <p). Let us mark 
the polar coordinates of the spin node as 9i, <f>i] from 
this set of coordinates, we can produce a 2J 7 '-dimensional 
complex vector: 

{Y lm {e h ^)}ff 1 . 

It is easy to see that if we apply a rotational symmetry 
operator g of the spinor condensate on this vector, we 
have: 



) = E4- is)Y hn (e jt <f> j ) 



(75) 

12 

where A\j (g) is the permutation operator from Eq. 
(74). The right-hand side of this equation indicates the 
rearrangement of the spin- nodes due to the symmetry op- 
erator. On the other hand, the left-hand side comes from 
our knowledge of the transformation rules for spherical 
harmonics, under rotations: namely, I, the total angu- 
lar momentum is invariant, and the different azimuthal 
angular momentum components mix under the transfor- 
mation. 

To connect the spherical harmonics with the z- 
representation, we need to construct a vector that will 
also transform with the phase e lXij< - 9 \ This requires 
that in addition to evaluating the spherical harmonics 
at the points (9i,<fii), we need to account for the phase 



factor when constructing the derived bases in the z- 
representation. This can be achieved by the following 
notion: instead of looking at the value of Y; m (#,</>), let 
us look at its derivative, which in our gauge convention 
can be written as 



BY, 



Ira \ 



dz* 



d_ 



._L d_ 

sin 9 dd 



Y lm (9, 0). (76) 



The denominator of the partial derivative dz* can be 
thought of as a small deviation, Sz*, from the mean- 
field spin node; it obeys the complex conjugate of the 
transformation rule in Eq. (74), so its inverse transforms 
in the correct way, using the phase e lXi ^ 9 \ Therefore, 
we finally have the connection between the symmetry of 
the spinor condensate and the representations of 50(3): 



(i) , x dY lm ,(6i, <j>j) 



0: 



dYu 



i) 



(77) 

What we achieved by making this connection is a way 
of constructing for each I > 1 (I = gives identically 
zero) partially reduced (albeit still reducible) represen- 
tations of the symmetry group at hand in terms of the 
z-paramctrization of small deviations from equilibrium. 
Let us denote the vectors we construct from Y\ m as: 



dY lm (0j, 

dz* J 



2F 



(78) 



These vectors are the simplest building blocks for the 
vibrational and rotational eigenvectors. 

As an example, consider the I = 1 states (p-states) ob- 
tained for the cyclic state of spin-two condensates. With 
the orientation for the cyclic state given in Sec. IV A, we 
obtain for 1 = 1, m = 0: 



111,0 oc {e i+ • x 3 }j =1 oc (1, 1, 1, 1) 



(79) 



which is the displacement vector for x 3 = z-axis rota- 
tion, as in Eq. (63). For m = ±1, as in atomic-orbital 
physics, it is useful to construct the p x and p y combina- 
tions, which are p x<y oc Yi i =F ii — i. For p x we obtain: 

Ui,i - ui _i cx {e i+ ■ £i}f =1 cx {1, -1, -1, 1}, (80) 

which is the rotation about the x\ = x-axis (up to an 
overall complex coefficient). In the same fashion we find 
that the p y combination is 

z M + zi _i cx {e l+ ■ x 2 }f =1 oc {1, -1, 1, -1} (81) 

which is corresponds to rotation about the x 2 = y-axis. 
To obtain the last mode, which is the optical vibration 
mode shown in Fig. 1, all we need is to find the vector 
of u which is orthogonal to the above three. 

The above analogy with atomic p-orbitals is not ac- 
cidental. Since we mapped vibrational modes to spher- 
ical harmonics, we also mapped the z-represcntation of 
spinor-condensate fluctuations to the Im representation 
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of atomic orbitals. In atomic orbital theory, we know 
that in the absence of rotational-symmetry breaking all 
m-states within the same I arc degenerate. But in the 
presence of a crystal field, this degeneracy is lifted. The 
effect of crystal fields on angular-momentum multiplets 
is very well-documented (see, e.g., [32]); we can now use 
this resource to directly find the eigenmodes of the spinor 
condensates. 

Let us demonstrate this principal again using the cyclic 
state. We have already shown that the I = 1 vibration 
modes correspond to rotations. Let us now consider the 
I = 2 states. Under the effect of a tetrahedral crystal 
field the electronic states split as: 

M -> J dxy ' dxz > dyz (82) 

[ CZ Z 2, d x 2_ y 2 

Now we can map back these atomic states to spinor- 
condensate oscillation modes. Starting with d xy cx Y2,2 — 
^2,-2 we find 

u 2>2 - u 2 _ 2 oc {1, 1, 1, 1} (83) 

which corresponds to uniform rotation about the z-axis. 
Similarly d xz oc Y 2 ,\ - Y2-1 and d yz cx F 2 ,i + Y 2 -i 
correspond to rotations about the y and cc-axis respec- 
tively. The two remaining orbitals are d z 2 cx Y 2 ,o an d 
d x 2_ y 2 cx Y 2 .2 + Y 2 ._ 2 . Since there are only four inde- 
pendent vectors, z, d z 2 and d x 2_ y 2 translate to the same 
zzm-vector: 

u 2 ,o oc u 2 . 2 + u 2 .-2 oc {1, 1, -1, -1} (84) 
which is exactly the optical mode shown in Fig. 1. 

V. CONCLUSIONS 

In this work, we applied the hydrodynamic descrip- 
tion developed in Ref. [1] to study the low lying exci- 
tations of the spinor condensate in the vicinity of the 
mean field-ground state. The dynamics of spinor con- 
densates close to the mean-field ground state is where 
their hidden point-group symmetry becomes most ap- 
parent and accessible. Using the spin-node formalism, 
and the parametrization of the spin-nodes in terms of a 
stereographic projection, we reduced the problem of find- 
ing the 2F spin- wave eigenmodes to a simple question of 
decomposing a representation of the appropriate point 
symmetry group to its irreducible representations. We 
also provided a simple recipe that allows the direct ex- 
traction of the condensate's spin- wave eigenmodes using 
the derivatives of the spherical harmonics, coupled with 
the knowledge of atomic orbital degeneracy lifting under 
a crystal field. Also, quite generally, we showed that all 
non-degenerate eigenvectors of the matrix Aij defined in 
Eq. (49) correspond to the eigenmodes of spinor con- 
densates, independent of the interaction (so long as it 
preserves SU(2) symmetry). 



More than any specific result, the current paper and 
Ref. [1] derive a new formalism to address high-spin many 
body systems. It is our impression that, by far, we have 
not yet explored all possible applications of this formal- 
ism. A simple example is the calculation of the spin- wave 
eigenmodes and energies of a spinor condensate which is 
locally at its ground state, but with its spin-nodes struc- 
ture rotated as a function of space. This can be done by 
combining the linearization of Sec. Ill with the general 
hydrodynamic description derived in Sec. II C. Similarly, 
our method of expanding about a mean-field ground state 
in terms of the z-variables could be readily applied to 
computing the leading instabilities in quantum-quench 
experiments (as in, for instance, Ref. [6] where spin-one 
quantum quench experiments were performed). The lin- 
earized Lagrangian derived in Sec. Ill C applies near any 
extremum of the spin interaction energy, Vi n t, even an 
unstable one. This can then be used to investigate the 
dynamics for short time-scales after a quantum quench. 

Another possible direction focuses on the form of the 
spin interaction energy T^„ t (ni, n 2 , . . . n 2 _F). In terms of 
the spin-nodes, the spin interaction energy must be a per- 
mutation symmetric function of the spin nodes. But the 
number of permutation symmetric scalars constructed of 
the spin nodes is limited. All such scalars must be 
constructed from tensors of the form: 

2F 
i=l 

where au = x,y, z is the space direction. Examples are: 

2F 2F 

n.^ • Xlj Tli,aL^ij,a. 

^=\ F (86) 

and so forth. This structure of the spin interaction may 
be used to construct generic phcnomcnological theories 
for spinor condensates and other high-spin many-body 
systems, along the lines of the construction of Landau 
free energy. 

The most interesting applications of the spin-node for- 
malism may arise when considering non-condensed spinor 
systems. Lattice insulators, both fermionic and bosonic, 
could also be parametrized using spin-nodes, and should 
exhibit magnetic mean-field states with hidden point- 
group symmetry as well. Similarly, we intend to consider 
spinor Fermi liquids using this formalism; such systems 
may have interesting magnetic instabilities into states 
with the same hidden symmetries as those arising in 
spinor condensates. 

In the challenging field of many body quantum sys- 
tems, often a new technical perspective on a problem 
may simplify it dramatically. In this paper we developed 
a formalism that seeks to do exactly that to the dynam- 
ics of spinor condensates - a topic of much current ex- 
perimental as well as theoretical interest. Our analysis 
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provides an economical representation, which allows for 
a direct, general, and easy calculation of many dynamic 
collective properties of spinor condensates. In addition, 
we hope that the developments presented here could be 
used in other challenging problems involving interacting 
quantum systems with high spin. 



This relation leads to a useful identity for the projec- 
tion operator 



\(Qi) 2F ) (T,n| 



(A4) 
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APPENDIX A: THE SPINOR BASIS {|T;fi)} AND 
THE MATRIX A 



In this Appendix, we will develop derive identities used 
for the projection operator V = 1 — \x) (x\ an d the sym- 
metry matrix A. Consider a particular spinor 



in) = \ci 1 n 2 ...n 2 F) 



(Al) 



where none of the spin nodes arc degenerate. Then from 
this we can construct a set of 2F states where one of the 
elements of |f2) is time-reversed {|T,;f2)}. Furthermore, 
we construct the set of 2F coherent states which are or- 
thogonal to |n) which are j (f2-) 2F ^)j. We note that 
these two sets of states satisfy reciprocal relations: 



<T 4 n| (o* 



2F, 



= A,; 5, 



(A2) 



where 



(2F)!n<%|0*) 



This relation can be immediately proved by expand- 
ing any state acting on the right in a basis of states 
| (p,f) 2F \ X, and any state acting on the left in a basis 

of states {|Tjf2)} (both which, in addition to the state 
|J1), form a complete basis of spinor states when the spin 
nodes are non-degenerate). 

Using these states, we will now proceed to derive an ex- 
pression for the inverse of the matrix Ay = (Tjf2| V \Tjtt) 
which exists when none of the spin nodes are degen- 
erate. We define B to be the matrix of the overlap of 
time-reversed coherent states (which will be shown to be 
the inverse of A) 



Ba 



m 2F (n 



,2F 



(A5) 



Consider the product of these matrices 



Y, B l3 A 3k = Y * ton| v |r fc n» 



(A6) 

We can then use the identity in Eq. (A4) to collapse the 
sum over j. This leads to 



£M>- W , W -fe (AT) 



(A3) and the proof is complete. 



[1] R. Barnett, D. Podolsky, and G. Refael, arxiv:0812.3403. 
[2] J. Stenger, S. Inouye, D. M. Stamper-Kurn, H. J. Mies- 

ner, A. P. Chikkatur, and W. Ketterle, Nature 396, 345 

(1998). 

[3] H. Schmaljohann, M. Erhard, J. Kronjager, M. Kottke, 
S. van Staa, L. Cacciapuoti, J. J. Arlt, K. Bongs, and 
K. Sengstock, Phys. Rev. Lett. 92, 040402 (2004). 

[4] M. S. Chang, C. D. Hamley, M. D. Barrett, J. A. Sauer, 
K. M. Fortier, W. Zhang, L. You, and M. S. Chapman, 
Phys. Rev. Lett. 92, 140403 (2004). 

[5] A. Griesmaier, J. Werner, S. Hensler, J. Stuhler, and 
T. Pfau, Phys. Rev. Lett. 94, 160401 (2005). 

[6] L. E. Sadler, J. M. Higbie, S. R. Leslie, M. Vengalattore, 



and D. M. Stamper-Kurn, Nature 443, 312 (2006). 
[7] M. Vengalattore, S. R. Leslie, J. Guzman, and D. M. 

Stamper-Kurn, Phys. Rev. Lett. 100, 170403 (2008). 
[8] S. Powell and S. Sachdev, Phys. Rev. A 76, 033612 

(2007). 

[9] R. Kanamoto, L. Carr, and M. Ueda, Phys. Rev. Lett. 
100, 060401 (2008). 

[10] Y. Kawaguchi, M. Nitta, and M. Ueda, Phys. Rev. Lett. 
100, 180403 (2008). 

[11] M. Taillefumier, E. Dahl, A. Brataasand, and W. Hof- 
stetter, arxiv:0901.1969. 

[12] A. Widera, F. Gerbier, S. Foelling, T. Gericke, O. Man- 
del, and I. Bloch, Phys. Rev. Lett. 95, 190405 (2005). 



13 



[13] F. Gerbier, A. Widera, S. Foiling, O. Mandel, and 
I. Bloch, Phys. Rev. A 73, 041602(R) (2006). 

[14] S. R. Leslie, J. Guzman, M. Vengalattore, J. D. Sau, 
M. L. Cohen, and D. M. Stamper-Kurn, arxiv:0806.1553. 

[15] M. Vengalattore, J. Guzman, S. Leslie, F. Serwane, and 

D. M. Stamper-Kurn, arxiv:0901.3800. 

[16] H. Saito, Y. Kawaguchi, and M. Ueda, Phys. Rev. A 76, 

043613 (2007). 
[17] A. Lamacraft, Phys. Rev. Lett. 98, 160404 (2007). 
[18] M. Uhlmann, R. Schiitzhold, and U. R. Fischer, Phys. 

Rev. Lett. 99, 120407 (2007). 
[19] B. Damski and W. H. Zurek, Phys. Rev. Lett. 99, 130402 

(2007). 

[20] S. Mukerjee, C. Xu, and J. E. Moore, Phys. Rev. B 76, 
104519 (2007). 

[21] G. I. Mias, N. R. Cooper, and S. M. Girvin, Phys. Rev. 

A 77, 023616 (2008). 
[22] R. W. Cherng, V. Gritsev, D. M. Stamper-Kurn, and 

E. Dernier, Phys. Rev. Lett. 100, 180404 (2008). 
[23] T. L. Ho, Phys. Rev. Lett. 81, 742 (1998). 

[24] T. Ohmi and K. Machida, J. Phys. Soc. Japan 67, 1822 



(1998). 

[25] C. V. Ciobanu, S. K. Yip, and T. L. Ho, Phys. Rev. A 

61, 033607 (2000). 
[26] M. Ueda and M. Koashi, Phys. Rev. A 65, 063602 (2002). 
[27] R. Diener and T.-L. Ho, Phys. Rev. Lett. 96, 190405 

(2006). 

[28] L. Santos and T. Pfau, Phys. Rev. Lett. 96, 190404 
(2006). 

[29] R. Barnett, A. Turner, and E. Demler, Phys. Rev. Lett. 

97, 180412 (2006). 
[30] S.-K. Yip, Phys. Rev. A 75, 023625 (2007). 
[31] H. Makela and K.-A. Suominen, Phys. Rev. Lett. 99, 

190408 (2007). 

[32] F. A. Cotton, Chemical Applications of Group Theory 
(Wiley, 1990). 

[33] G. Herzberg, Molecular Spectra and Molecular Structure, 
Volume II (Van Norstrand Company, 1945). 

[34] C. J. Pethick and H. Smith, Bose-Einstein Condensation 
in Dilute Gases (Cambridge University Press, 2002). 



